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ABSTRACT 

We report on the surprisingly high metalhcity measured in two absorption systems at 
high redshift, detected in the Very Large Telescope spectrum of the afterglow of the 
gamma-ray burst GRB 090323. The two systems, at redshift z — 3.5673 and z = 3.5774 
(separation Av w 660 km s~^), are dominated by the neutral gas in the interstellar 
medium of the parent galaxies. From the singly ionized zinc and sulfur, we estimate 
oversolar metallicities of [Zn/H] = +0.29 ±0.10 and [S/H] = +0.67 ±0.34, in the blue 
and red absorber, respectively. These are the highest metallicities ever measured in 
galaxies at z > 3. We propose that the two systems trace two galaxies in the process of 
merging, whose star formation and metallicity are heightened by the interaction. This 
enhanced star formation might also have triggered the birth of the GRB progenitor. As 
typically seen in star-forming galaxies, the fine-structure absorption Si ii* is detected, 
both in GO and Gl. From the rest-frame UV emission in the GRB location, we derive 
a relatively high, not corrected for dust extinction, star-formation rate SFR « 6 Mq 
yr~^. These properties suggest a possible connection between some high-redshift GRB 
host galaxies and high-z massive sub-millimeter galaxies, which are characterized by 
disturbed morphologies and high metallicities. Our result provides additional evidence 
that the dispersion in the chemical enrichment of the Universe at high redshift is 
substantial, with the existence of very metal rich galaxies less than two billion years 
after the Big Bang. 

Key v^rords: Cosmology: observations — galaxies: ISM — gamma-ray burst: indi- 
vidual (GRB 090323). 



1 INTRODUCTION 

The chemical enrichment history of galaxies is one of the 
key diagnostics of the evolutionary tracks of our Universe. 
The main element production factories since the Big Bang 
are stars, which transform light elements into heavier ones. 
Understanding the cosmic chemical evolution requires de- 
tailed spectroscopic observations and theoretical models. For 

* Partly based on observations collected at the European South- 
ern Observatory under ESO proposal No. 082.A-0693. 
t E-mail: savaglioOmpe. mpg.de 



a long time, the high-redshift Universe was explored by an- 
alyzing absorption systems in galaxies distributed along the 
sight lines to high-redshift QSOs. Absorbers with high H I 
content (A'hi ^ 10^*^ cm~^), called damped Lyman-alpha 
(DLA) systems (Prochaska & Wolfe 1997), have predom- 
inantly neutral gas (A'h — A^hi, temperature T < 1,000 
K). Their metallicity is generally low at high redshift, way 
below 1/10 solar at z ~ 4, and approaches solar in the 
nearby Universe (Wolfe, Gawiser, & Prochaska 2005; Per- 
oux et al. 2003). More recently, it has been possible to in- 
vestigate the chemical enrichment in galaxies using emis- 
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Figure 1. The optical VLT/FORS2 spectrum of the afterglow of GRB 090323. Indicated are the main absorption lines associated with 
the two systems GO and Gl at ^ 3.57. Many other lines, associated with intervening systems are present, but not marked. Znll 
absorption is identified in Gl. The rod lino at the bottom is the spectrum of the Earth atmosphere. The two stars indicate the telluric 
absorption. 



sion lines, whose fluxes represent the integrated signal emit- 
ted by hot gas {T ~ 10, 000 K) in H ll regions of star 
forming galaxies. About solar metallicity was detected in 
2.1 < z < 2.4 galaxies (Shapley et al. 2004). Similarly, the 
mass-metallicity relation detected in Lyman-break galaxies 
(LBGs) at 3.0 < z < 3.7 indicates that the most massive 
objects (M* ~ 10" M0) can be metal rich (Maiolino et al. 
2008; Mannucci et al. 2009). The emerging picture is that 
the spread in metallicity is large at any redshift, certainly 
much larger than what one would expect based on DLAs in 
QSO spectra alone. 

DLAs are regularly detected in spectra of gamma-ray 
burst (GRB) afterglows. Long GRBs are the most lumi- 
nous explosions associated with the death of massive stars, 
and, as regularly done with QSOs, are used as background 
bright point-like sources for the investigation of the inter- 
stellar medium (ISM) in high-2 galaxies. The main differ- 
ences to QSOs are that GRB afterglows can be brighter, 
are short lived and happen in star-forming regions. QSO- 
DLAs trace primarily metal-polluted galaxy haloes, where 
most of the volume of a galaxy is located and where dust 
extinction is less important (Fynbo et al. 2008). The first 
chemical analysis of absorption lines in the DLA detected in 
a 2; ~ 2 GRB afterglow spectrum revealed a high (about so- 
lar) metallicity (Savaglio et al. 2003). This neutral-ISM dom- 
inated GRB-DLA in a high-z galaxy, the GRB host, opened 
a new view our understanding of the chemical enrichment 
of the Universe. In particular, the about solar abundance 
called into question the idea that low metal pollution in 
galaxies is the rule at high redshift. The typical GRB-DLA 



metallicity is indeed larger than in QSO-DLAs (Savaglio et 
al. 2003; Vreeswijk et al. 2004; Berger et al. 2006; Prochaska 
et al. 2007b; Fynbo et al. 2008), leading to the conclusion 
that high-z GRB hosts show properties more similar to UV 
bright star forming, LBG like, galaxies, than QSO-DLAs. 



Here we report the discovery of two DLAs at 
z = 3.5774±0.0005 (henceforth called GO) and z = 
3.5673±0.0003 (Gl), separated in the velocity space by 
Ri 660 km s~^, revealed in the rest-frame UV afterglow spec- 
trum of GRB 090323 (Chornock et al. 2009; Cenko et al. 
2011; McBreen et al. 2010). The two absorbers (probing the 
Universe when it was 13% of its present age) could belong 
to a single galaxy, to two distinct galaxies in close proximity 
and perhaps in the process of merging, or to two physically 
distant and unrelated galaxies. In the latter case, if the red- 
shift difference is attributed to cosmic expansion alone, it 
would correspond to a separation of ~ 5.8 Mpc in comoving 
coordinates. 



In this paper we analyze the spectrum of the 
GRB 090323 afterglow to probe the chemical enrichment of 
GO and Gl. Observations are described in §2, the spectral 
analysis and heavy-element measurements in §3 and §4, re- 
spectively. The final discussion and conclusions are in §5 
and §6, respectively. Throughout the paper we adopt an 
Ho — 70 km s~^, Qm ~ 0.3, SIa ~ 0.7 cosmology (Spergel 
et al. 2003). 
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Figure 2. The metal absorption lines associated with the neutral 
gas in GO and Gl, at Av = +660 km s~^ and A?) = km s~^, 
respectively. The dotted line is the noise spectrum. 

2 OBSERVATIONS 

GRB 090323 was detected and localized by the Gamma- 
ray Burst Monitor (GBM) and the Large Area Telescope 
(LAT) onboard NASA's Fermi Gamma-ray Space Telescope 
(Michelson et al. 2008) on 23 March UT (Ohno et al. 
2009). Observations by Swift's X-ray Telescope (XRT; Ken- 
nea et al. 2009) lead to the discovery of the afterglow and 
enabled longer-wavelength follow-up by the Gamma-Ray 
burst Optical/Near-infrared Detector (GROND; Updike et 
al. 2009; McBreen et al. 2010), Palomar 60-inch Telescope 
(Cenko et al. 2009), and Gemini-S (Chornock et al. 2009). 

Optical spectroscopy of the afterglow started on 2009 
March 24.25 UT (1.2 d post-burst) using the FOcal Reducer 
and low dispersion Spectrograph 2 (FORS2) mounted at 
the 8.2 m ESO-VLT Ul telescope in Paranal, Chile, when 
the AB 7?-band magnitude of the afterglow was « 18.6 
(McBreen et al. 2010). We obtained 5 x 900 s integrations us- 
ing the 600RI grism (5300 A- 8450 A coverage) and 6 x 900 s 
integrations using the 600z grism (7500 A- 10,000 A cover- 
age) with a long slit of I'.'O width. The bluer sequence started 
at 06:05 UT and ended at 08:13 UT, and was alternated with 
the redder sequence which started at 06:53 UT and ended 
at 09:01 UT the same night. The average seeing was ~ 0'.'60 
resulting in effective spectral resolutions of ~ 4.6 A FWHM 
at 6300 A, and 7» 4.8 A FWHM at 8100 A, or « 210 km s"^ 
and « 150kms~^, respectively. 

The data were reduced with standard IRAF routines, 
and spectra were extracted using an optimal (variance- 
weighted) method. Spectro-photometric flux calibration 
was carried out against observations of the standard star 
LTT6248 obtained in the same night. Corrections for slit 
losses, due to finite slit width and for Galactic foreground ex- 
tinction, were applied. The individual exposures were com- 



bined and the final spectrum is shown in Figure [T] In order 
to search for spectral time variability (in particular in the 
fine-structure transition Sill*), we also merged the first and 
second halves of the exposures individually and compared 
the strengths of the absorption lines. No variability was de- 
tected. 



3 COLUMN DENSITIES OF ABSORPTION 
LINES 

Details of some metal lines associated with the low- 
ionization species in Gl and GO are shown in Figure (2] GO 
is the reddest absorption system detected in the spectrum 
of GRB 090323. However, given the possible strong interac- 
tion between the two, we cannot be sure that GO is asso- 
ciated with the galaxy hosting GRB 090323. Column densi- 
ties (Tables [T] and [2)) are estimated by using the apparent 
optical depth (AOD) method (Savage & Sembach 1991) and 
curve of growth (COG) analysis (Spitzer 1978). These meth- 
ods are routinely applied when high-resolution and signal- 
to-noise spectroscopy (for which the best fit of Voigt profiles 
is possible) is not available. The COG analysis (Figure [3]) 
is robust for weak lines, when the linear approximation can 
be applied (e.g., rest-frame equivalent width Wr < 0.2 A 
and effective Doppler parameter 6 > 30 km s~^; Spitzer 
1978). For stronger lines, the COG gives reliable results if 
several lines with different oscillator strengths of the same 
ion are detected. This idea is supported, for instance, by the 
spectral analysis of the afterglow of GRB 081008 by D'Elia 
et al. (2011), who have shown that column densities derived 
from a low resolution spectrum with the COG are consistent 
with those derived from Voigt profile fitting of a high res- 
olution spectrum. The AOD method is used for lines with 
moderate to low equivalent widths (EWs), provided that 
the effective Doppler parameter is not very low {h > 10 
km s~^; Prochaska et al. 2007a). COG values (adopted for 
our final abundance determination) are generally consistent 
with AOD ones (Tables [T] and . The AOD gives at times 
slightly lower values, indicating that residual saturation is 
not always taken into account. Prochaska (2006) has shown 
that column densities of metal lines derived with the COG 
can be underestimated. If this problem affected our mea- 
surements, the true metallicities would be even larger than 
those reported here. 

The H I column densities are derived from the Lyman-a 
spectral absorptions (Figure |4]). These are strongly blended 
and show the characteristic damping wing typical for col- 
umn densities > 10^^ cm~^. The Voigt-profile best fit gives 

iVHI = 1020'72±0.09 ^^-2 ^ ^Ql9.62±0.33 ^^^-2 

Gl and GO, respectively (at fixed redshifts, as derived for 
the metal systems). As indicated by Figure |4l there is virtu- 
ally no lower bound to A'hi in GO. However, the presence of 
strong absorption lines associated with heavy-element low- 
ionization transitions is an indication that A^hi cannot be 
very low. The upper limit for the metallicity in Gl is set by 
the requirement that here the metallicity cannot be higher 
than 10 times solar. 

In Figure [5] we also display the velocity profiles of high- 
ionization lines associated with GO and Gl (e.g., Nv, Siiv, 
and Civ). Together with GO and Gl, we mark a possible 
weak absorption system ~ 500 km s~^ blueward of Gl, with 
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Figure 3. Curve of growth analysis of low-ionization species in Gl (Figure[3ji) and GO (Figure[3l)). In the x-axis label, f\ indicates the 
oscillator strength of the atomic transition. The best fit and Icr uncertainty are the continuous and dotted curves, respectively, while 
the straight line shows the linear approximation. The column density of the medium-ionization species Al III for both systems is also 
included. For Gl, the Doppler parameter 6 ~ 84 km s~^ has been derived from the best fit of Sin, which is consistent with the best-fit 
Doppler parameter of Sill* and Aim. For all other ions, this is assumed to be 84 km . Similarly for GO, the Doppler parameter b ~ 90 
km s^^ has been derived from the best fit of Sill. The best fit curve of growth for Aim gives a very uncertain b ~ 48 km s~^, but the 
best-fit column density is not very different if we force b = 90 km s^^, as for Sill. 



the tentative identification of NvA1239 and ClvA1548 ab- 
sorptions. The VLT spectra also reveal a wealth of strong 
metal absorption lines associated with absorbers in galax- 
ies distributed along the GRB sight line, at redshift z > 1. 
These are at z ^ 1.697, 2.101, 2.448, and 3.379, and not 
considered further in this study. 



4 HEAVY-ELEMENT ABUNDANCES 

In DLAs, the ionization correction can generally be ne- 
glected (Dessauges-Zavadsky et al. 2003; Jenkins 2009), 
and ions with ionization potentials just above the ioniza- 
tion energy of hydrogen dominate all other states (e.g., 
iVpoii ~ iVpe, A^znii ~ A^zn, A^s II ~ A^s)- The chemical 
state of the gas is thus derived using the simple approxi- 
mation [X/H] = log(A''xii/A''Hi) - log(A'^x/A''H)0 (Pettini et 
al. 1997). Among the many detected lines, the interesting 
ZnnA2026 line in Gl at z = 3.5673 is detected with high 
confidence (Figure [2| thanks to the high signal-to-noise in 
the spectrum. 

The Hi column density of GO is relatively low (A''hi = 
2^019.62 f-Qj. ^ijicij the ionization correction could be 

necessary. The ionization correction is generally smaller than 
0.2 dex in DLAs with 19.5 < log(iVHi/cm"2) < 20.3, and 
smaller for larger Hi columns (Dessauges-Zavadsky et al. 
2003; Peroux et al. 2007). The low ionization in our sys- 
tems is supported by large column densities of low-ionization 
metal lines: Kanekar et al. (2009) found an anti-correlation 
between the gas temperature and metallicity of DLAs. The 
relatively low column density of Aim (Table [2]) additionally 



indicates low ionization (Meiring et al. 2008). Although we 
can only infer an upper limit for Al ll, the expected value can 
be estimated from other ions with similar ionization poten- 
tial, from which we conclude that A'aiii/A'aiiii > 10. Some 
heavy elements (e.g., iron or nickel) can be partly locked up 
in dust grains, and measured column densities (which are 
for the gas phase only) can be very different from the total 
column density. Sulfur and zinc are mainly in the gas phase, 
so that uncertain dust depletion corrections are less critical. 

Figure [B] shows the abundances of Si and S relative to 
Zn as a function of Znll column density for GO and Gl, 
other GRB-DLAs and QSO-DLAs. For Gl we derive [S/Zn] 
= -0.33 ±0.04 and [Si/Zn] = -0.72 ±0.06. Anomalous rela- 
tive abundances have been found before in high- 2: absorption 
line systems. In QSO-DLAs, significant [S/Zn] < is ob- 
served for about half of all systems where these two elements 
are detected. Cooke et al. (2011) report on a QSO-DLA with 
[S/C] < —1.51. The metal column densities, derived from 
a composite spectrum of 60 GRB afterglow spectra, indi- 
cate [S/Zn] = -0.94 (Christensen et al. 2011). Nissen et al. 
(2004) report low S-to-Zn ratios in QSO-DLAs with respect 
to Galactic stars, from which a different star-formation his- 
tory is hypothesized. Nissen et al. (2007) were able to lower 
[S/Zn] values in Galactic stars by using non-LTE modeling. 
Vladilo et al. (2011) suggest that the observed deficit of S 
is indicative of a larger dust depletion with respect to Zn in 
dust rich sight lines. These anomalies are not totally under- 
stood, another contributing factor might be a theoretically 
very uncertain element yields in early generations of massive 
and low-metallicity stars. 

The ZnllA2026 line in Gl is contaminated by Mgl and 
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Table 1. Element abundances of Gl at z = 3.5673 (Nyii = 1020 72±0.09 cm-2). 
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"An estimated small contamination by MglA2026 {Wr — 0.031 A) is subtracted from the measured EW. 
CrllA2026 is very weak and can be neglected (see text). 
''Zn IIA2062 is not detected because of blending with strong telluric lines. 
'^Contaminated by SllA1250 of GO. 
Contaminated by OlA1302 of GO. 



Table 2. Element abundances of GO at z = 3.5774 (A^hi = lois ^^^icaa ^m'^). 
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An upper limit to Zn IIA2062 is not estimated because of blending with strong telluric lines. 
''Contaminated by S IIA1253 of Gl. 



Cm. CrllA2026 is weak and can be neglected. MglA2026 is 
rarely detected at high spectral resolution and columns are 
in the interval A^Mgi = ]^q12.3-13.2 ^^-2 -p^ estimate the 
MglA2026 contamination, we considered EWs measured in 
the GRB afterglow composite spectrum of Christensen et 
al. (2011). The EWs of metal lines in Gl and the compos- 
ite are similar, in particular for weaker lines {Wr < 1.1 A). 
From this, we assumed a possible Wr ~ 0.8 + 0.3 A for 
MgIA2852 in Gl (this is not covered by our spectrum). The 
error is given by the dispersion of the EWs for Gl and GRB 



composite around the equality relation (80% of the lines). 
Thus, from MgIA2852 and the COG (we assumed fo = 84 
km s~^, as for other lines), we estimated Wr = 0.0311q qj4 
A for MgIA2026, which we subtracted from the ZnllA2026 
EW (Table [1)). The CrllA2056 line is in a noisy part of the 
spectrum (sky emission contamination). We set a conserva- 
tive limit Afcrii < 10^*" cm"^, which gives [Cr/Zn] < -0.3, 
explainable by a possible depletion of Cr into dust grains. 

A Wolf-Rayet (WR) star origin was advocated to ex- 
plain the presence of double absorption features in GRB 
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Figure 4. The two panels show the FORS2 spectrum of 
GRB 090323 around the Lya region. The two systems at 2 = 
3.5673 (Gl) and z = 3.5774 (GO) are indicated by the two verti- 
cal tick marks. Upper panel: the dashed-blue and red-dotted lines 
indicate the individual Lya profiles with AfHi,Gl = lO^"'^'^^'''''® 
cm~^ and A'^hi.GO = \Qi19.62±0.33 cm~^, respectively. From these, 
we derive [Zn/H] = -1-0.29 ± 0.10 and [S/H] = -1-0.67 ± 0.34. 
Lower panel: the red solid line is the total best-fit profile. The 
two dotted lines represent two distinct solutions contained in 
the 90% confidence interval. These are A^hIjCi 



10^ 



and A^HI.GO = 10^ 
iVHl.G0 = 10" *« c 
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spectra. This could also explain the presence of the two 
strong nearby absorbers GO and Gl in GRB 090323. The 
GRB would take place within the dense, fast-moving shell 
of low-ionized gas. This scenario is disfavored (Chen et al. 
2007), because the GRB strongly ionizes the circumburst 
environment out to a kpc distance scale (Vreeswijk et al. 
2007; D'Elia et al. 2009; Ledoux et al. 2009), and thus large 
column densities of low ionized gas in the wind, as indicated 
by our observations, are unlikely. Local WR stars indeed 
do not show strong low-ionization absorption lines. Further- 
more, even if hydrogen remained neutral in the blast, and 
20 Mq of hydrogen were ejected spherically to a distance of 
~ 10 pc, we would expect a Hi column density of 6 x 10^* 
cm~^, which is well below the measured value. 

The association of GO and Gl with a single star forming 
region and its wind is also unlikely. We considered a standard 
SN remnant scenario and, under the assumption of a shell 




-4000 aOOO 2000 -'1000 -3000 2000 
velocity (km s"') velocity {km s">) 

Figure 5. High-ionization absorptions in GRB 090323. The two 
strongest absorbers are marked by the two vertical dashed lines. A 
possible third weak absorber, shown by ClvA1548 and NvA1239 
features, at z = 3.5600 (Ad ~ 500 km s~^ blueward of Gl) is also 
marked by the vertical dotted line. The horizontal dotted line is 
the noise spectrum. 



thickness of ~ 10% of the shell radius, compute the swept- 
up mass. We consider ISM densities over a wide range, from 
1 cm~^ to 1500 cm~"^, and compute the shell radius for the 
time when the shell velocity differential reaches the observed 
660 km s~^. The resulting swept-up mass is of the order of 
1000 M© . A single SN cannot enrich this amount of swept-up 
ISM mass to a metallicity a few times solar. Even a situation 
of triggered star formation, with dozens SNe during a small 
time interval cannot change this picture, since each SN has 
to sweep up its 800-1000 M© to slow down to the observed 
velocity spread of 660 km s~^. 

4.1 Fine-structure absorption in GO and Gl 

Absorption lines associated with excited and metastable 
transitions are regularly reported in GRB afterglow spectra 
(Ledoux et al. 2009). The first excited transitions J3/2 asso- 
ciated with the fine-structure species Si 11* are detected both 
in GO and Gl (Tables [1] and and Figure [2]) with column 
densities A^siii. = lO^^^ i^io o^ ^^^-2 ^q13.82±o.05 ^^-2^ 

respectively. No other excited or metastable transitions are 
observed (some lines are likely blended with strong ground- 
level lines). In QSO-DLAs, Cll*, but not Sill*, is regularly 
observed (Wolfe et al. 2003). Figure [7] (left panel) shows 
the column densities of several ions (including Sill*) in Gl 
and the star-forming galaxy cB58 at z = 2.72 (Pettini et al. 
2002). The similarity between the column densities in Gl 
and cB58 suggests that the physical conditions of the neutral 
gas in Gl are not dissimilar from those of high- 2; star-forming 
galaxies. In addition, Pettini et al. (2002) measured a rel- 
atively high metal enrichment: [S/H] = -0.36 ± 0.09. The 
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Figure 6. Silicon-to-zinc and sulfur-to-zinc relative abundances 
(upper and lower panel, respectively) as a function of Zn II column 
density. Limits for GO (with N'Znll < 10^^'^ cm~^) and values for 
Gl (with Afznil = lOl^ *^ cm-2) are the blue filled circles. Other 
GRB-DLAs are represented as red squares, QSO-DLAs are black 
dots. Dashed lines indicate solar relative abundances. 



Detected absorption lines indicate that the contribution 
to extinction of GO can be neglected (column densities are 
weaker than in Gl). To estimate the possible extinction in 
Gl, we use the method described in Savaglio & Fall (2004). 
As a reference value, we take the rate of extinction per col- 
umn of neutral gas derived for the Milky Way (Bohlin et al. 
1978): 

iVui 



< MW 



0.53 X 



102icm~ 



(1) 



This choice of Milky Way extinction rate is supported by a 
relatively low value of the iron-to-zinc ratio [Fe/Zn] = — 1.5 
in Gl, not too dissimilar to the typical ratio in the ISM of 
the Milky Way (Savage & Sembach 1991). Eq.[T]can also be 
expressed in terms of the column of metals A*""'^* , by assum- 
ing that the rate of extinction in the optical per column of 
dust Av/Af''""* is proportional to Av/iV""°' (Savaglio & Fall 
2004), or: 



<^ = 0.53 



jymct 

ivrmct 
^^0,21 



(2) 



where Nq°2i is the metal column in a gas cloud with A^'h i 



0.53 for 



10 cm~ and solar metallicity. Eq.[2]gives Ai 
a gas cloud with zinc column density A'zn ~ 10^ 
Our Znll column density for Gl is A'znii = io^^-ST gjjj-a^ 
which means that Ay ~ 0.5 is expected. This optical ex- 
tinction is higher than that measured from the SED {Ay = 
0.1 — 0.2), and the typical Ay derived from GRB afterglow 
SEDs (Zafar et al. 2011). The relatively flat GRB 090323 af- 
terglow SED cannot be explained by a possible large column 
of dust (see also a similar conclusion in Savaglio & Fall 2004) , 
unless a sufficiently grey extinction is assumed (Stratta et 
al. 2004). However, a grey extinction is not confirmed by the 
mean dust extinction law derived from a sample of 17 GRB 
afterglow sight lines with X-ray-to-optical SEDs (Schady et 
al. 2011). 



right panel of Figure [7| shows the comparison between the 
ionized gas in Gl and close to the QSO APM 08279-1-5255 
at z = 3.9 (Srianand & Petitjean 2000). The Sill*-to-Sill 
ratio in GO and Gl is small (0.005 and 0.01, respectively), 
much lower than in APM 08279-1-5255, where this is ~ 0.3, 
and indicating a high ionization of the gas near the mas- 
sive black hole. The typical Sin*-to-Sin ratio measured in 
GRB-DLAs is > 0.01. Following the modeling of Srianand 
& Petitjean (2000), we conclude that the absorbing gas in 
GO, Gl and cB58 has a low density (< 100 cm~^) and/or a 
low temperature (T < 10^ K). 

4.2 Dust depletion and extinction 

McBreen et al. (2010) and Cenko et al. (2011) derived the 
optical dust extinction in the sight line of GRB 090323 by 
using the SED of the afterglow. This is relatively fiat red- 
ward of the Lya, indicating a moderate Ay = 0.1 — 0.2. We 
alternatively estimated Ay from the metal column densi- 
ties of several elements with different refractory properties. 
These give information on the depletion of metals onto dust 
grains and thus the expected dust extinction. 



5 DISCUSSION 

In Gl we found the most distant Znll absorption ever de- 
tected up to now, being followed by a DLA at 2 = 3.39 in 
the spectrum of QSO 0000-2620 (Molaro et al. 2000). Zinc 
is widely considered to be the most reliable indicator of the 
metal abundance of DLAs in high-redshift galaxies (Pettini 
et al. 1997), and preferred over more abundant elements (e.g. 
iron, silicon or oxygen) because it is much less depleted to 
the dust component of the ISM, and because the spectral 
analysis is less affected by saturation. The column density 
iVznii = 10^^-^''='='' °* cm"^ in Gl is also among the highest 
ever inferred, with other systems of comparable density all 
residing at 2 < 2.9. The large column densities in Gl would 
indicate the galaxy to be massive and/or metal rich. At this 
high redshift, this would make Gl a rare object. 

In Figure |8] we show the redshift distribution of metal- 
licities as derived in GRB-DLAfO and QSO-DLAs. DLAs 
detected in QSO spectra are distributed between the QSO 



^ For more on GRB-DLA metallicities, see Berger et al. (2006), 
Fynbo et al. (2006), Savaglio (2006), Prochaska et al. (2007b) and 
Rau et al. (2010). 
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Figure 7. The comparison between the column densities of Gl and cB58 (left panel), and Gl and QSO APM08279+5255 (right panel). 
The straight lines on the left and right panels indicate the x = y relations. This shows the similarities between the two DLAs in Gl and 
cB58 and the strong difference between Gl and narrow absorption lines close to QSO regions. Column densities for cB58 are recalculated 
by Savaglio (2006) and consistent with those reported by Pettini et al. (2002). 



and the observer, and not directly related to the QSO it- 
self. On the other hand, GRB-DLAs are generally associated 
with the galaxy hosting the burst. The average metallicity 
and dispersion in 17 GRB-DLAs in z = 2.0 — 4.5 (excluding 
GO and Gl) is log(Z/Z0) = -0.98 ± 0.64, whereas this is 
log(^/Z0) = -0.83 ± 0.76 if we include GO and Gl. The 
average metallicity in 186 QSO-DLAs in the same redshift 
interval is log(.Z/Zo) = -1.39 ± 0.61 (Figure [8|. Zinc and 
sulfur abundances in Gl and GO are [Zn/H] = -1-0.29 ± 0.10 
and [S/H] = -1-0.67 ± 0.34, respectively (blue filled circles 
in Figure [8]). Uncertainties are dominated by errors on Hi 
columnfl TVhi = lO^o-^^io.og ^^-2 ^ ^q19.62±o.33 

cm~^ for Gl and GO, respectively. However, these (and 
thus the metallicities) are not independent of each other. 
A higher A^hi in one system would imply a lower A'hi in 
the second. If we consider one absorber for the Lya, we de- 
rive Nni = 102°-^''±° °* cm"^ and z = 3.5690 (the redshift 
is not fixed), which is equivalent to the total HI columns 
of GO plus Gl (A^Hi = 102° '^5±° °^ cm-2). From this, we 
get a total metaUicity [Zn/H] = +0.25 ± 0.09 and [S/H] 
= -1-0.07 ± 0.08. Regardless, both absorbing systems ex- 
hibit super-solar metallicities, well above QSO-DLAs and 
all other GRB-DLAs (Figure [8]). The high abundances in 
GO and Gl argue against an explanation invoking a natural 
metal enhancement in star-forming environments associated 
with the birth sites of GRB progenitors (Prochaska et al. 
2007b; Fynbo et al. 2008): it would be very unlikely for the 
two metal-rich regions, widely separated in velocity space, to 
be located both nearby GRB 090323. Adopting a very con- 
servative approach, we estimate a lower limit to the metallic- 
ity by assuming that the total HI column density is in either 



^ The upper bound to A^hi of GO is artificially set to give a 
metallicity 10 times the solar value. 



of the two systems. These would give [Zn/H]> -fO.25 ± 0.09 
and [S/H]> -0.47 ± 0.09. 

The highest metallicities in the local Universe are gener- 
ally found in massive quiescent elliptical galaxies, for which 
further metallicity enrichment is prevented by the lack of 
cold gas, necessary for star formation. High metallicities 
are reached very rapidly (Matteucci 1994), which means, 
given the old age of the stellar population, at high redshift. 
The mass-metallicity relation in galaxies suggest that the 
typically low stellar mass of z < 2 GRB hosts (of the or- 
der of ~ 10^'^ M0) is not in contradiction with the low 
metallicities measured (Savaglio et al. 2009; Levesque et 
al. 2010). What happens at higher redshift is not totally 
understood. The high metallicity in GO and Gl suggests 
high stellar masses. Model predictions in Figure [8] indicate 
M. > 10" M© (Savaglio et al. 2005). Although these are 
based on measurements from emission lines, which indicate 
the average metallicity of the entire galaxy, these must be 
correlated with those measured in a DLA along one sight 
line, with some dispersion. The stellar component where GO 
and Gl have origin was detected with the multi-band imager 
GROND (Greiner et al. 2008) and the measured AB mag- 
nitudes are r' = 24.87 ± 0.15, i' = 24.25 ± 0.18, J > 22.7, 
H > 22.2 and K > 21.6 (McBreen et al. 2010). Unfortu- 
nately, the detection in the optical only is not constrain- 
ing the stellar mass because it probes the rest-frame UV, 
where the mass-to-light ratio can vary by a factor > 100. 
The emission (not resolved) has a size < 6 kpc, but this 
is not constraining either, because galaxies (both early and 
late types) in the distant Universe are several times smaller 
than local galaxies with similar masses (Buitrago et al. 2008; 
van Dokkum et al. 2008; Damjanov et al. 2009; Toft et al. 
2009; Mosleh et al. 2011). From the optical detection, we 
can estimate the total (integrated over GO and Gl) star- 
formation rate (SFR) by assuming that the rest-frame UV 
emission originates in the hot gas where stars form, and us- 
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ing the relation between the rest-frame UV flux (for instance 
at A = 1500 A) and SFR in GRB host galaxies (Savaglio et 
al. 2009). From the observed fluxes at z ~ 3.57, we estimate 
Lisoo = 4 X 10^° erg s'^ A-\ and a non-dust corrected SFR 
= 6.4 Mq yr~^. This value is relatively high, if one considers 
that the (unconstrained) dust attenuation in the galaxy can 
easily boost the SFR. 

The GO and Gl absorbers could be in a single mas- 
sive galaxy, with high metallicity and velocity dispersion. 
They could also be two galaxies tidally disrupted early on, 
in a metal-rich environment, for instance in a galaxy clus- 
ter, whose metallicity in the intracluster medium can rise 
above the solar value (Sanders & Fabian 2006). If the gas 
mass in GO and Gl is not dominating over the total baryonic 
galaxy mass, further residual star formation would increase 
the metallicity over time by a small quantity. The system 
would eventually evolve into an elliptical galaxy after a ma- 
jor merger event, without dramatically affecting our view 
of the chemical state of the local Universe, where over-solar 
metallicities are possible in quiescent M, > 10^'' galax- 
ies. Alternatively, observational biases may have prevented 
us from detecting objects like GO or Gl in larger numbers. 
Star formation proceeded very rapidly in the past, which 
would have led to a very rapid increase in the metallicity 
and in the exhaustion of gas. What would remain today 
would be very metal-rich galaxies escaping observations be- 
cause star formation would have ceased early on, and ten 
billion years later the system would contain mainly old and 
faint stars. 

Double absorbers are not rare in GRB afterglow spec- 
tra. We statistically studied the fraction of these strong dou- 
ble absorption systems with small physical separation in 
GRB afterglow spectra (Savaglio et al. 2002; Klose et al. 
2004; Ferrero et al. 2009; Page et al. 2009) and found 5 out 
of 40 spectra (~ 12%). For comparison, we studied a sample 
of about 500 QSO-DLAs, and found 18 (about 4% of the 
total) close pairs (or triplets) at a distance in velocity space 
< 10,000 km s"\ The excess in GRB-DLAs suggests that 
an enhanced star formation due to interaction and merger of 
galaxies could be a key mechanism triggering high-z GRB 
events (Lopez-Sanchez & Esteban 2009; Chen 2011), and 
not, as commonly assumed, the low metallicity of the pro- 
genitor (MacFadyen & Woosley 1999). At z < 1, Mannucci 
et al. (2011) suggested that the low metallicity measured 
in GRB hosts can be explained by the fundamental mass- 
metallicity relation: galaxies with similar stellar mass have 
higher star formation rate if more metal-poor. Campisi et 
al. (2011) studied this problem using cosmological simula- 
tions of GRB hosts and also concluded that a high SFR 
is likely the fundamental parameter correlated with GRBs. 
The task now is to seek similar indications at high redshift. 
High metallicities and interactions are suggested by another 
class of high— z galaxie^: the rare dusty and exceptionally 
luminous sub-mm galaxies (SMGs; Blain et al. 2002). Metal- 
licities are estimated by converting the sub-millimeter emis- 

^ Solar and super-solar metallicities at 2 > 3.5 were measured 
before, in the proximity of QSOs, by using narrow absorption and 
broad emission line systems (Srianand & Petitjean 2000; Dietrich 
et al. 2003; Jiang et al. 2007). However, these systems contain 
1/10 or less of the total baryonic content of their host galaxies 
and have a marginal effect on the global metal content. 
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Figure 8. Metallicities as a function of redshift of cold ISM 
in galaxies. Gl and GO metallicities in GRB 090323 (left and 
right blue dots, respectively, artificially separated by Az = 0.06 
for clarity) are given by [Zn/H] = -1-0.29 ± 0.10 and [S/H] = 
-1-0.67 ± 0.34, respectively. Metallicities of other GRB-DLAs and 
QSO-DLAs are red squares and black crosses, respectively. GRB- 
DLA metallicities are collected from the literature (Prochaska et 
al. 2007a; Savaglio et al. 2009), except for GRB 080210 (z = 2.641 
and logZ/Z0 = -0.93; De Cia et al. 2011), GRB 081008 (z = 
1.968 and logZ/Zg = -0.52; D'Elia et al. 2011), GRB 090926 
{z = 2.106 and logZ/Zg = -1.89; Rau et al. 2011), and 
GRB100219A (z = 4.667 and logZ/Zg = -1.10; Thone et al. 
2011). The dashed line is the linear correlation for QSO-DLA 
points. Solid curves are metallicities expected for star-forming 
galaxies with different stellar masses, derived by using the empir- 
ical mass-metallicity relation and its redshift evolution (Savaglio 
et al. 2005), after proper normalization at z = 0. Metallicities are 
derived using Znll, Sin, Sll, and Fen for GRB-DLAs, and Znn, 
Sin and Fen for QSO-DLAs, after dust-depletion correction. The 
upper bound for GO is set by assuming that its metallicity can- 
not be higher than 10 times solar. Some GRB-DLAs do not have 
error-bars due to lack of reported uncertainties on the measured 
column densities. 

sion (which traces molecules and dust) into mass of metals 
(Blain et al. 2002), for instance assuming, as found in the 
local Universe, that the dust mass is about 40% of the to- 
tal mass of metals (Franceschini et al. 2000). Metallicities 
derived from the nebular emission lines in SCUBA galaxies 
at 1.4 < z < 2.7 are slightly below solar (Swinbank et al. 
2004). A connection between some high-z GRB host galaxies 
and SMGs is supported by the high luminosities {L > 10^^ 
Lq) found in 20% of the GRB hosts observed in the ra- 
dio with SCUBA (Berger et al. 2003). From the optical-to- 
radio SED of 0.8 < z < 1.5 GRB hosts, Michalowski et 
al. (2008) suggest that these galaxies can be similar to the 
hotter (thus less luminous) sub-sample of SMGs. In GRB 
hosts, the dust-hidden (not detected in UV/optical) star 
formation is revealed in the sub-mm, from which a more 
likely SFR > 100 Mq yr"'^ is measured (Castro Ceron et al. 
2006; Michalowski et al. 2008; Hunt et al. 2011). This can 



10 S. Savaglio et al. 



be a factor 10 to 100, or more, higher than measured in the 
UV/optical (Savaglio et al. 2009). If this applies also to GO 
and Gl, for which the UV emission indicates SFR > 6 M© 
yr~^, a SFR higher than 100 Mq yr~^ is possible. The role 
of dust is still far from been sufHciently understood, but a 
better insight is now provided by early and deep SED obser- 
vations of complete samples of GRB afterglows (Greiner et 
al. 2010; Kriihler et al. 2011). It is found that the majority 
of those defined as dark GRBs were instead dimmed by dust 
along the sight lines. 

The situation might be different in the nearby Universe. 
In a recent work, Yates et al. (2011) report that mergers 
of massive galaxies at low redshift are much more frequent 
when the metallicity is low, in contrast with what is pro- 
posed here for the z ~ 3.6 Universe. The chemical evolution 
of interacting galaxies has been investigated by Perez et al. 
(2011) and Torrey et al. (2011) with cosmological simula- 
tions. A firm conclusion is that the star formation history 
in massive mergers is complex and different from isolated 
galajcies. According to Torrey et al. (2011), gas rich systems 
have higher metallicity after the interaction than metal poor 
ones. 



6 SUMMARY 

We report the analysis of the heavy-element enrichment of 
two DLAs at z = 3.5673 and z = 3.5774 detected in the 
VLT/FORS2 spectrum of the afterglow of GRB 090323. One 
or both absorbers, separated in velocity space by about 660 
km s~^, are in the vicinity of the GRB event itself. Which 
one is closer is not clear. If the velocity difference is purely 
due to the cosmological expansion, the comoving distance 
between the two absorbers would be 5.8 Mpc. The blue and 
the red absorption systems (galaxies Gl and GO, respec- 
tively) show very strong metal lines, from which we derive 
zinc and sulfur over-solar abundances in the ISM: [Zn/H] 
= +0.29 ± 0.10 and [S/H] = +0.67 ± 0.34, respectively. The 
two systems could be physically separated, but because of 
the high metal abundance in both, we propose the possi- 
bility that the two systems are interacting and in the pro- 
cess of merging. Ground optical detections have not enough 
resolution to resolve the components. Nevertheless, the UV 
rest-frame image gives a high star formation rate SFR > 6 
M0 yr~^, not corrected for dust extinction. The double ab- 
sorbers, the high metallicity at high z and high SFR suggest 
that the two galaxies are related to massive sub-mm galax- 
ies. If this is the case, the interaction might be one important 
condition through which high-z GRB events occur. The sit- 
uation is different from what is seen at 2: < 2, where GRB 
hosts generally show relatively low metallicities and stellar 
masses. Considering that the global star formation rate of 
the Universe moves from big to small systems, going from 
high to low redshift, our discovery supports the idea that 
massive star formation is the main mechanism triggering 
GRBs and not metallicity. 

The discovery of metal rich galaxies at z > 3 is impor- 
tant for our understanding of the cosmic structure formation 
and heavy-element enrichment from a metal free era to the 
present time abundances. GRBs are certainly special events 
that might occur in special environments. However, through 
follow-up spectroscopy of afterglows, we can highlight re- 



gions of the Universe that are hard to detect in other ways. 
Most of the gas probed by our observations is perhaps in a 
galaxy nearby the GRB. If galaxies as those discovered by 
us are common, an unrecognized population of metal-rich 
galaxies might exist, which could significantly alter our un- 
derstanding of the cosmic chemical evolution. The challenge 
now is to find such remnants. 
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